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1. Introduction 

Lattice Quantum Chromo Dynamics (QCD), denned on a space-time lattice, has been 
our best and most reliable tool to extract non-perturbative physics of the strongly 
interacting theory for over a decade now. In many cases, such as the hadron spectra or 
the weak decay constants, the focus of lattice QCD has now shifted to obtaining more 
precise results. In the fields of strong interest to this conference, namely Relativistic 
Heavy Ion Collisions and the possible transition to the new Quark-Gluon Plasma (QGP) 
phase, also there are similar examples : the quark-hadron transition temperature pQ, 
T c , the Wroblewski Parameter |2j, A s and the the equation of state pQ. Already a lot 
of reliable theoretical information on these physical observables of interest has been 
predicted by lattice QCD and again the focus is now for better precision. Due to the 
strong influence of dynamical fermions on thermodynamics and the nature of the phase 
transition, one still has some way to go before precise prediction for realistic world of two 
light (u and d) and one heavy (s) flavour become available. Nevertheless, the continuing 
progress in lattice QCD makes it perhaps a reachable goal. I shall be unable to review 
all the important recent lattice results in the time available. I have therefore chosen to 
focus on those with thrust on something new, either conceptually or quantitatively. In 
particular, I shall cover the new results on the T-/i# phase diagram and the physically 
interesting question of J/^-dissolution/persistence. An interesting theoretical issue 
these days has been the comparison of predictions of conformally invariant theories 
with lattice QCD results. I shall show one such comparison along with new results for 
speed of sound. 

2. QCD Phase Diagram 

The first principles exploration of the QCD phase diagram has the lattice QCD partition 
function as its focal point. Assuming three flavours of quarks, u, d, and s and denoting 
by fif the corresponding chemical potentials, it is defined by, 



Note that the quark mass rrif and the corresponding chemical potential /// enter only 
through the Dirac matrix M for each flavour. The baryon and isospin chemical potentials 
are related to those for quarks by fiB = fi u + + Us and ^ = fi u — \id- 

Progress in unravelling the phase diagram using lattice QCD has been slow due to 
so-called "fermion problem": for nonzero fi, the determinant in eq. is complex. This 
renders both analytical and numerical investigations very difficult. Recently, progress 
was made by a key observation. Following the phase transition from the /i# = axis into 
the hb~T plane is relatively easier. Various methods have been employed to achieve this. 
In addition to the original re-weighting method [3] in two couplings, imaginary chemical 
potential jlj and Taylor series expansion have been used to find and locate the critical 
end point, expected for a world with two light flavours on the basis of a variety of 




f=u,d,s 



(1) 
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model considerations [U] . I shall describe below our results obtained by using the Taylor 
series method. One of its primary advantages is the ease of taking the continuum and 
thermodynamic limit which is mandatory for all such lattice computations. This is 
especially difficult for the re-weighting method [3] due to the exp(— AS) factor it has 
which becomes exponentially small for bigger lattices. Furthermore, the discretization 
errors propagate in an unknown manner in a re-weighting computation, making reliable 
estimates in continuum limit very difficult. A much better control of the systematic 
errors is, on the other hand, achieved in the method we employ. We study volume 
dependence of the Taylor coefficients at several temperature T to i) bracket the critical 
region and then to ii) track its change as a function of volume. As shown in Figure 1, 
one expects to pin down the critical region much more precisely this way in going from 
a smaller volume V\ to a bigger one V2. 

T 




Figure 1. Volume dependence of the critical region. 



Various quark, baryon and isospin densities and the corresponding susceptibilities can 
be obtained from eq. ((H) as: 

_Td\nZ Td 2 \nZ 

Higher order susceptibilities, Xn u ,n d , can be similarly defined as n u {rid) partial derivatives 
of liaZ with respect to yU n (/^)- Using them, the pressure P has the following expansion 
in fi: 

AP _ P(fM,T) P(0,T) _ 1 / //.. \ n « 1 / // j\ n d 



n~<4. r 7~'4 r 7~'4 / j A,*iu,*id 

11 1 n u ,n d n i 



(3) 



Quark number susceptibilities in (j2J) are crucial for many quark-gluon plasma 
signatures which are based on fluctuations in globally conserved quantities such as 
baryon number or electric charge. Theoretically, they serve as an important independent 
check on the methods and/or models which aim to explain the large deviations of the 
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lattice results for pressure P(/i=0) from the corresponding perturbative expansion. This 
has been discussed elsewhere in details. Here we will be concerned with using them 
with the higher order terms to learn about the critical end point. Indeed, a series for 
baryonic susceptibility xb can be constructed using them, whose radius of convergence 
gives the nearest critical point 7 . Limited by the finite (and even small) number of 
terms available, one uses standard tests for series convergence on them to do the best 
one can. Successive estimates for the radius of convergence can be obtained from these 
using r n+2 = \Jxb/Xb +2 - We use terms up to 8th order in /x, i.e., the estimates from 
2/4, 4/6 and 6/8 terms in eq.([3*)). It may be noted here that a similar construction of a 
series goes through for the off-diagonal susceptibility, xn as well- The ratio Xn/Xio can 
be shown to yield the ratio of widths of the measure in the imaginary and real directions 
at n = 0, i.e., a measure of the difficulty of the 'fermion sign' problem mentioned above. 
With some care, this argument can be generalized to nonzero \i. 

2.1. Our Results 




Figure 2. The radius of convergence as a function of N s (m w N s ). The left panel 
shows results obtained by comparing the 4th and 6th order at T/T c = 0.9 (pentagons), 
0.95 (circles) and 1 (triangles). The right panel shows results obtained at T/T c = 0.95 
using the ratio of the 4th and 6th order NLS (boxes) and the 6th and 8th order NLS 
(circles) . 

Our results jZj were obtained by simulating full QCD with two flavours of dynamical 
staggered quarks of mass m/T c = 0.1, where T c is the transition temperature for zero 
chemical potential, on lattices of size 4 x iVj*, with N s = 8, 10, 12, 16 and 24. R- 
algorithm [Hj with a trajectory length of one unit of molecular dynamics time on N s = 8 
was used for simulations of full QCD. The trajectory length was scaled with N s on larger 
spatial lattices. Earlier investigations have pinned down the physical parameters for our 
simulations at zero chemical potential [9, : m p /T c = 5.4 ±0.2 and m^/m p = 0.31 ±0.01. 
We used the canonical methods to fix the lattice spacing a in physical units for which 
we made runs on symmetric 16 4 lattices. Their details as well as other details about 
evaluation of the higher order susceptibilities can be found in Here I shall be 
brief and state only that our simulations were made at T/T c = 0.75(2), 0.80(2), 0.85(1), 
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0.90(1), 0.95(1), 0.975(10), 1.00(1), 1.05(1), 1.25(1), 1.65(6) and 2.15 (10) with a typical 
statistics of 50-100 in the units of (maximum) autocorrelation length. 

Using terms of up to 8th order in \i in the expansion in eq. (JHJ), we obtained 
estimates for the radius of convergence shown in Figure 2. The left panel shows r 6 as 
a function of N s at three temperatures whereas the right panel exhibits r$ and r§ at 
T/T c = 0.95. The latter also shows on top the spatial lattice size in units of inverse 
pion length. As seen from the left panel of the figure, our estimate of the radius of 
convergence for the smaller lattices with N s ~ 10 is consistent with the critical end 
point estimate of the re-weighting method jSj with appropriate size, i.e, N s m v = 3-4. 
Note that m^/irtp = 0.31 in that well. Figure 2 shows, however, large finite size 

effects for these small lattices and a strong change in the estimate around N s m n ~ 6. 
Such finite size effects are expected theoretically [TH] and have been seen in other cases 
such as chiral condensates [TT] and lattice QCD evaluations of structure functions [T2"] . 
A key consequence of our use of larger volumes than before is that we find a shift in 
the critical end point to //b/T ~ 1 — 2, as is evident from Figure 2. Recently, similar 
computations but up to 6th order have been reported on a large lattice with N s m n ~ 15, 
and were also discussed in this conference [IB]. These use, however, large quark masses, 
and consequently a a large m n /m p ~ 0.7. It will be interesting to see if lowering their 
quark mass leads to results similar to ours. 

2.2. Measure of Sign Problem 

0.1 , , , , 
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Figure 3. The ratio of quark number susceptibilities %2 and xn on 4 x 16 3 (circles) 
and 4 x 24 3 (boxes) lattices. 



Figure 3 displays the ratio of the two susceptibilities xn an d X20 as a function 
of temperature on two different lattice sizes. This ratio, which is a measure of the 
seriousness of sign problem, shows how the increase in both its value and fluctuations 
make the simulations progressively more difficult as one lowers the temperature to T c 
and below. The success of the current round of finite density lattice simulations can 
also be traced from this figure to the fact that all current efforts, regardless of the 
details of their methods, attempt to explore the phase diagram starting from T c , the 
transition point at zero chemical potential. An important check one can make in the 
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Taylor series expansion technique is that each coefficient in the Taylor expansion must 
be volume independent. Since the actual evaluation of these terms entails cancellation 
of divergences, this is indeed a nontrivial check which has been made successfully [Zj. 
In fact, care is needed in extracting any equation of state at finite chemical potential 
from these coefficients : a systematic study on different volumes (spatial lattices) with 
a subsequent extrapolation to infinite volume is necessary to ensure that some peak-like 
behaviour in, e.g., X40; is rea l an d does not give spurious results for thermodynamic 
quantities. 

3. Speed of Sound 
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Figure 4. Specific heat C v and 4e (left) and the speed of sound c s as a function of 
temperature for quenched QCD. See text for more details. 



An important physical property of the strongly interacting matter at finite 
temperature, which has been computed [H] since long using lattice techniques, is the the 
speed of sound. Its relevance to many aspects of heavy ion physics, such as the elliptic 
flow or the hydrodynamical studies can hardly be overemphasized. A related quantity of 
interest is the specific heat which has been related to event-by-event fluctuations of the 
transverse momentum (pr)- There are various measures of these fluctuations currently 
in use in the literature. It may actually be better to relate them to the specific heat in 
order to make them both more physical as well as uniform in convention. 

A recent development in this direction on the lattice is a new way to obtain these 
by relating them to the temperature derivative of an anomaly measure A/e = (e — 3P)/e 
[To] . Using lattices with 8, 10, and 12 temporal sites (38 3 x 12 and 38 4 lattices) and 
with statistics of 0.5-1 million iterations, e, P, s, C 2 S and C v were obtained in continuum 
using this differential method at 2T C and 3T C for quenched QCD. Sizable deviations from 
the corresponding ideal gas results were found although C v was consistent with 4e, as 
expected in a weak coupling scenario. On the other hand the result for entropy density 
at 3T C also seemed to agree with a strong coupling (super Yang-Mills) prediction jTHj : 
s/ So = f(g 2 N c ), where f(x) = § + § C(3)(2a;- 3 / 2 ) + ■ • ■ and s = |vr 2 A^ c 2 T 3 . However, 
curiously it too failed at 2T C . 
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New results for the specific heat and the speed of sound, obtained by using a 
new method, called the t-method, were presented at this conference in a poster |17j . 
While the differential method ^H] for obtaining thermodynamic quantities has lesser 
systematic errors and is convenient to use, it was abandoned in favour of the integral 
method [19] due to the former's problem of yielding negative Pressure. It was shown 
that this problem of negative pressure can be completely avoided even in the differential 
method by writing the derivatives with respect to temperature and volume using a 
parameter t : the old (Bielefeld) method [IHj corresponds to t = 1 and leads to negative 
pressure but t = results in positive pressure always. Figure 4 shows these new results 
of [TJj for the specific heat (the left panel) and the speed of sound (right panel). While 
the corresponding pressure is not shown here, it is positive even around T c , an input 
which went in getting the speed of sound (squares), shown in Figure 4. Note that these 
continuum results have been obtained by simulating lattices with many lattice spacings 
and then extrapolating to zero unlike the results for the integral method (dotted line) 
in Figure 4 which are for N t = 8. The interesting and noteworthy features are : i) 
the specific heat does tend to peak near T c and is widely different there from the weak 
coupling expectation of 4e, ii) for T > 2T C , this weak coupling expectation is fulfilled 
but the resulting value is very different from that for an ideal gas, iii) speed of sound 
is about 10 % smaller than the ideal gas for higher temperatures but drops close to 
T c , and iv) the results are in good agreement with the integral method but could differ 
significantly in the neighbourhood of T c . Further investigations are needed to see if c 2 s 
dips close to zero near T c . 

4. Persistence of J/ip 




0246246 135135 
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Figure 5. Spectral functions [2Dl of the vector J/ip and scalar \ mesons as a function 
of energy at various temperatures. See text for more details. 

An intriguing result on the J/ijj and similar heavy quarkonia has recently emerged 
PIH I2T] from the lattice QCD investigations. As is well known, Matsui and Satz [22] 
proposed J/ip suppression as a signal of quark-gluon plasma. Quarkonium potential 
model calculations with an ansatz for temperature dependence of the heavy quark- 
antiquark potential lead to the result that J/ip and Xc states dissolve in QGP by 
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1.1T C . The impressive results from the CERN NA50 [23| did suggest an anomalous 
suppression of J/ip in Pb-Pb collisions, bringing forth detailed theoretical questions 
on this suppression mechanism and the above mentioned ansatz. In particular, one 
approach taken by the lattice experts [2H1 EE] has been to study the spectral function of 
these states in QGP to look for their dissolution from first principles. The recognition 
of the maximum entropy method (MEM) as a reliable tool to extract the spectral 
functions from temporal lattice correlators at zero temperature lead to its application 
at finite temperature. Caution is, however, called for in doing so since nonzero and high 
temperatures are obtained by making the temporal lattice smaller whereas MEM works 
better for longer temporal lattices. 

Figure 4 shows the results of the Bielefeld group 20J, obtained by simulating lattices 
from 48 3 x 12 to 64 3 x 24 in sizes in the quenched approximation of no dynamical quarks. 
The vertical bars are drawn to indicate the size of the error on the integrated spectral 
function in the area shown by the horizontal bars. One observes that the Xc does seem 
to indeed dissolve by 1.1T C , as expected from naive potential model considerations. 
However, the J/ip and rj c persist up to 2.25 T c and are gone only at 3T C . Similar results 
have been obtained by other groups |2I| as well, although there seems to be some 
disagreement on the precise value of the temperature at which the J/ip peak vanishes. 

While these results have to be still confirmed in a full QCD simulation with light 
dynamical quarks, they already can have strong consequences for the J/ip signal. Recall 
that only about 30-40 % of the J/ ip come through \ decays ( with a substantial fraction 
from X2 which has not yet been investigated ), thus the predicted suppression will be 
rather small until temperatures close to 3T C are reached (or energy densities 81 times 
that of critical energy density is reached), which seems to be a tall order for the current 
round of CERN and BNL experiments. An interesting observable consequence could, 
however, be a change of the suppression pattern as a function of y/s or the energy 
density reached. 

5. Summary 

The QCD phase diagram in T — plane has begun to emerge. Since all results so 
far have been obtained on N t = 4, i.e, very coarse, lattices, several quantitative issues 
related to continuum limit still have to wait. It may be noted that the existence of 
equally legitimate [24. different prescriptions of introduction of the chemical potential 
on the lattice can, and do, lead to even larger lattice artifacts on such coarse lattices. 
Nevertheless, it is very heartening to find that different methods lead to a similar 
qualitative picture. 

Our results, using Taylor expansion, permit us to take the thermodynamic limit, 
unlike other methods. Our investigation of volume dependence suggests N s m n > 6 
approximates the thermodynamic volume limit reasonably. The [Ib/T of the critical 
end point, identified from the radius of convergence using up to, and including the 
8th order terms, drops strongly around that volume. We find that yUg/T ~ 1 — 2 is 
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indicated for the critical point with a corresponding temperature of ~ 0.95T C , where T c 
is the transition temperature for the baryonfree case. 

1.1 , , , , , , 



1 - 




Figure 6. The QCD phase diagram with a freeze-out curve superimposed. The filled 
circle denotes the estimate of the critical point which has been obtained in [7|. The 
open circle is an earlier estimate from [2] using smaller lattices and nearly the same 
quark mass. 

Figure 6 shows our results, together with the earlier determination j^j on small 
volumes but similar quark mass, i.e, the pion mass. A freeze-out curve obtained 
by treating a resonance gas as ideal, has been superimposed on the Figure. From the 
marked values of the CM energy, a/s, per nucleon needed to reach that point on the 
curve, one notes that the critical end-point may be in the reach of future RHIC energy 
scans. 

I also presented new continuum results on the speed of sound and the specific heat 
for quenched QCD. A new method, which revives the earlier differential method but 
without its drawback of negative Pressure was used to obtain these results. Interesting 
quantitative agreement was found with strong coupling predictions from super Yang- 
Mills theory at 3T C but not at 2T C . 

Finally, the intriguing persistence of J/ip in QGP was discussed. The MEM 
technique has been used to obtain the spectral function at finite temperature. While it 
showed the peak in the x spectral function to go away by 1.1T C , the peak in the vector 
spectral function persists up to 2.25-3T c . Extension of these simulations to incorporate 
dynamical quarks would be crucial to be sure of the phenomenological consequences of 
these results. 
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